""" Implements zero temperature molecular dynamics with a simple integrator. """


import numpy


class NVE:

    def __init__(self, speck, potential, dt = 0.1):
        self.speck = speck
        self.potential = potential
        self.dt = dt
        self.speck.v = speck.r * 0.0
        
    def step(self):
        self.potential.force(self.speck)
        self.speck.v += (self.speck.f / self.speck.m) * self.dt
        self.speck.r += self.speck.v * self.dt
    
